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ABSTRACT 

We introduce the Anharmonic Oscillator Symmetry Model to describe vibrational ex- 
citations in molecular systems exhibiting high degree of symmetry. A systematic procedure 
is proposed to establish the relation between the algebraic and configuration space formu- 
lations, leading to new interactions in the algebraic model. This approach incorporates the 
full power of group theoretical techniques and provides reliable spectroscopic predictions. 
We illustrate the method for the case of Da/j-triatomic molecules. 
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1. Introduction 



Spectroscopic techniques represent one of the most important tools in modern chemical 
analysis. 1 ^ In particular, the molecular vibrational degrees of freedom are studied by means 
of Infrared and Raman spectroscopy. 2 ) It is necessary, however, to rely on theoretical 
models in order to interpret the data, which in turn refines the models in a feedback cycle. 

The study of molecular vibrational excitations is carried out by taking into account 
different degrees of approximations and theoretical assumptions. The simplest way to 
study the molecular energy spectra is by means of a Dunham expansion. 2 ) This method, 
however, does not provide wave functions, and consequently does not allow the calculation 
of physical properties such as transition intensities. On the other hand there are ab initio 
calculations, where an exact solution of the Schrodinger equation is attempted. In practice, 
the molecular Hamiltonian is usually parametrized as a function of internal coordinates 
and the potential modeled in terms of force field constants, 3 ) which are determined either 
through calculations involving the molecular electronic states for several configurations 4 ) 
or empirically, by the fitting of experimental data. 5 ) While for diatomic and triatomic 
molecules very accurate information on force field constants is available, 6 ) this is not the 
case for polyatomic molecules, due to the large size of their configuration spaces. It is 
thus important to develope alternative calculational methods in order to describe complex 
molecules for which ab initio calculations are not feasible. 

Algebraic models attempt to provide such alternative techniques. In 1981 an algebraic 
approach was proposed to describe the roto-vibrational structure of diatomic molecules, 7 ) 
subsequently extended to linear tri- and four- atomic molecules 8 ) and certain non-linear 
triatomic molecules. 9 ) Although these results were encouraging, the model could not in 
practice be extended to polyatomic molecules, for which it is necessary to incorporate 
the underlying discrete symmetries. This difficulty can be surmounted by treating the 
vibrational degrees of freedom separately from the rotations. In 1984 Van Roosmalen 
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et al proposed a U(2)-based model to describe the stretching vibrational modes in ABA 
molecules, 10 ) later extended to describe the stretching vibrations of polyatomic molecules 
such as octahedral and benzene like molecules. 11 -* Recently the bending modes have also 
been incorporated to the framework, which was then applied to describe C2„-triatomic 
molecules 12 ^ and the lower excitations of tetrahedral molecules, 13 ^ using a scheme which 
combines Lie-algebraic and point group methods. In a different approach, it has also been 
suggested to use a U(k + 1) model for the k = 3n — 3 rotational and vibrational degrees 
of freedom of a n-atomic molecule. This model has the advantage that it incorporates all 
rotations and vibrations and takes into account the relevant point group symmetry, 14 ^ but 
for larger molecules the number of possible interactions and the size of the Hamiltonian 
matrices increase very rapidly, making it impractical to apply. 

The algebraic formulations have no doubt proved useful, but several problems re- 
mained, the most important of which is the absence of a clear connection to configuration 
space traditional methods, which in turn makes their significance difficult to gauge. A 
related problem is the lack of a systematic procedure to construct the physically mean- 
ingful interactions in the algebraic space. In this paper we address both these issues and 
introduce a general model for the analysis of molecular vibrational spectra, which we call 
the Anharmonic Oscillator Symmetry Model (AOSM). We shall show that it is possible to 
construct algebraic operators with a well defined physical interpretation and in particular 
the interactions which are of special relevance for the description of the degenerate modes 
present in systems exhibiting high degree of symmetry. These are derived in a procedure 
that takes full advantage of the discrete symmetry of the molecule and that provides all 
possible terms in a systematic fashion. As a bonus, a clear-cut connection is established 
between the algebraic scheme and the traditional analyses based on internal coordinates, 
which correspond to the harmonic limit of the model. 

As a test for this approach we apply the AOSM to three "D 3 /i-triatomic molecular 
systems, namely Na^, Be3 and H3". We have chosen this set of molecules because they 
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cover a wide range of behaviors. Whereas Na^ is very harmonic, H3" displays a strongly 
anharmonic spectrum, while the Be3 cluster has an intermediate behavior. Since these 
small molecules can be very well described by means of ab initio calculations, 15 ' 16 -* we 
again emphasize the aim of this work. We shall establish an exact correspondence be- 
tween configuration space and algebraic interactions in the harmonic limit of the U(2) 
algebra. This general procedure not only allows to derive the interactions in the AOSM 
from interactions in configuration space, but can also be applied to cases for which no 
configuration space interactions are available. The Da/i-triatomic molecules constitute the 
simplest systems where degenerate modes appear and where the new interactions in the 
model become significant. The application of these techniques to more complex systems, 
such as tetrahedral molecules, is presented elsewhere. 17 ^ 

In the next section the structure of the model is presented. Section 3 is devoted to the 
construction of the symmetry adapted normal basis, which is the most appropriate one to 
carry out the diagonalizations. In Section 4 we propose a new set of interactions which 
have physical interest and suggest the need to construct operators associated to the E 
mode. In Section 5 we describe the general procedure to derive the algebraic interactions 
from those appearing in configuration space and in Section 6 we introduce the AOSM in 
order to derive all algebraic interactions from symmetry considerations. In Section 7 we 
apply the model to H3", Be3 and Na^, while in Section 8 we present our conclusions and 
discuss some future developments of the model. 
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2. Algebraic Model 



The model exploits the isomorphism of the U(2) Lie algebra and the one dimensional 
Morse oscillator 

ft 2 ^72 

" = -2?5? + fl ( e " T - 2e "'> • (21 > 

whose eigenstates £ can be associated with £7(2) D SO{2) states. 18 ) In order to see how 
this isomorphism comes about, consider the radial equation 

-^---r- + ^+r 2 J^(r) = (iV + l^(r) , (2.2) 

which corresponds to a two-dimensional harmonic oscillator (in units where h = p = e = 1) 
associated to a £7(2) symmetry algebra. 19 ^ By carrying out the transformation 

r 2 = (N + l)e- p , 



equation (2.2) transforms into 



d V^V 2 '-2e-) 



0(p) = -m^ip) , (2.3) 



dp 2 v 2 

which can be identified with (2.1) after defining x = pd and multiplying by h 2 /2pd 2 , 
provided that 

D = J^(N + l)> , (2.4 o) 
S = . (2.46) 

In the framework of the U{2) algebra, the operator N corresponds to the total number 
of bosons and is fixed by the potential shape according to (2.4), while m, the eigenvalue 
of the SO (2) generator J z , takes the values m = ±N/2, ±(N - 2)/2, .... The Morse 
spectrum is reproduced twice and consequently for these applications the m-values must 
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be restricted to be positive. In terms of the U(2) algebra, it is clear from (2.3-4) that the 
Morse Hamiltonian has the algebraic realization 

In addition, the U{2) algebra includes the raising and lowering operators J+, J_, which 
connect different energy states in (2.3), while the angular momentum operator is given by 
J 2 = ±N(N + 2), as will be shown below. 

The Morse Hamiltonian (2.5) can be rewritten in the more convenient form 

H = AH M = ^[(J + J_ + J_J+)-N] , (2.6) 

where we have used the relation J| = J 2 — |(J + J_ + J-J+) and added the constant 
term in order to place the ground state at zero energy. The eigenf unctions of the 
Hamiltonian (2.7) are associated to the U(2) D 0(2) chain, and are given by 

U(2) D SO (2) 

I I 
\[N] , v> 

with 



\[N],v>= yj (J-y\[N},0> (2.7) 

where N is the total number of bosons fixed by the potential shape (Eq. (2.4a)) and v 
corresponds to the number of quanta in the oscillator. Both, iV and v, are related with 
the usual labels j and m of the U(2) and 50(2) groups, by means of 19 ) 

N = 2j 

v = j — m . (2-8) 

The parameters iV and A appearing in (2.6) are related to the usual harmonic and 
anharmonic constants uj e and x e u e used in spectroscopy. 7 ) This is seen by substituting the 
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operator J z in (2.5) by its eigenvalue. In terms of v, the corresponding energy expression 
takes the form 

Em = —^(N + 1/2) + A(N + l)(v + 1/2) — A(v + 1/2) 2 , (2.9) 
from which we immediately obtain 

u e = A(N + l) , 

X e UJ e = A . (2.10) 

Thus, in a diatomic molecule the parameters A and N can be determined by the spectro- 
scopic constants u e and x e u e . 

We now consider a particular molecular system. We start by assigning a U l (2) algebra 
to each relevant interatomic interaction. 13 '' In figure 1 we show the U l (2) assignment for 
D3/,-triatoniic molecules. All relevant operators in the model are then expressed in terms 
of the generators of the molecular dynamical group, which is given by the product 

U\2) <g) U 2 (2) <g> U 3 (2) . (2.11) 

A simple realization for these generators can be given in terms of the number operator Ni 
and the operators J M; i 

{N{, Jx,ii Jy,i-i J z,i\i i = 1) 2, 3 

with 

N i = 8 \si + t\u (2.12a) 
J x ,i = 2 (A* + 4t), Jy,i = \{tU - sju), J Zji = -(s\ Si - t\u) , (2.126) 
where s\(si) and t\{tj) are bosonic operators satisfying the usual commutation relations 



7 



All other commutators vanish. The operators (2.12b) satisfy the usual angular momentum 
commutation relations. Computing Jf = + + J^, we find 




(2.13a) 



from which the identification 

ji = Ni/2 , (2.136) 

is readily made. One can show directly from (2.12) or from (2.13a) that [J M) i, Ni] = 0. The 
set (2.12b) thus defines the SUi(2) subalgebra of Ui(2). Since N = Ni will remain fixed in 
our applications, we shall sometimes refer to SUi(2) instead of Ui(2). The specific boson 
realization (2.12) was given for reasons of clarity, but is not necessary for the subsequent 
developments of the model. 

Formally, while the vibrational symmetry group of the X3 molecules is V^h, m practice 
it can be reduced to V s due to the in-plane restriction. Since we are assigning a number to 
each bond it is more convenient to work with the symmetric group S3, which is isomorphic 
to V 3 through the generator identification 

C 3 ^(123) , 
C, a ~(23) , 

as indicated in Fig. 1. The Hamiltonian of the system is then expanded in terms of the 
dynamical group generators (2.12), provided that we impose its invariance with respect to 
the symmetry group S3. In order to explain the main features of the algebraic model we 
start by considering a simple form for the Hamiltonian, restricted to two body interactions 
which preserve the total number of quanta V = £V Vi, where each Vi is defined as in (2.8), 

H = A^t 1 + f E H i3 + - 2 f: % , (2-14) 
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where the operators Hf 1 , and Vj are defined as 

Hf 1 = ^(j +i j_ i + j_ i j +i )-^ , (2.15a) 

Hij = 2JoiJoj 7^ ' (2.156) 

% = (J +l J-j + J-J+j) . (2.15c) 

The first term in (2.14) corresponds to three equivalent Morse oscillators (2.6), while the 
two terms and Vij correspond to interactions diagonal in the chains associated to the 
couplings 

SU {i) (2) <g> SU U) {2) D SO il) (2) <g> SO U) {2) D SO {lj) {2) , (2.16a) 

SU {l) (2) ® SU {j) (2) D SU&)(2) D 50^(2) , (2.166) 

respectively. The notation SU^i^) indicates the usual angular momentum coupling of 
the SU®(2) and SU^{2) states. 

The basis arising from three couplings of the form (2.16a) is refered to as the local 
basis, since the Morse oscillators are diagonal when the three SO 1 (2) algebras are well 
defined. 10 ' 19 ) It should be noted that for most calculations, higher order terms are required 
in the Hamiltonian (2.14) in order to attain higher accuracy, as we shall see in the following 
sections. The physical interpretation of these interactions will also be explained. 

Once we have established the form of the Hamiltonian, we need a basis to carry out 
its diagonalization. Since the Hamiltonian is invariant under the symmetry group S3, its 
eigenfunctions span irreducible representations (irreps) of this group for any given basis. 
It is convenient, however, to define a physical basis in order to classify the states with the 
usual normal mode labels, as well as to simplify the calculations. 
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3. Symmetry Adapted Normal Basis 



The simplest basis to diagonalize the Hamiltonian (2.14) is the one associated to the 
local mode chain 19 ) 

t/W(2) <g> C/( 2 )(2)(g) U^(2) D SOM(2)® SO( 2 \2)® SO^(2)d SO (2) 
I [iVi] , [N 2 ] , [N 3 ] ; Vl , v 2 , v 3 ; V > , 

where below each group we have indicated the eigenvalues that label their irreps. Explicitly 
this basis is given by, 

\[N 1 ],[N 2 ],[N 3 ];v 1 v 2 v 3 >=\[N 1 ];v 1 > \[N 2 ];v 2 > \[N 3 ];v 3 > (3.1) 

where [[ATj]; Vi] are given by (2.7). The index Vi corresponds to the number of quanta in the 
i-th oscillator, which is related to the eigenvalues ji and rrii of the J 2 and J Zji operators 
by means of expressions (2.8) and (2.13) 

Vi = ji - rrii . (3.2) 

Since 

ji = y , mi > , (3.3) 

we find 

u i = 0,l,2,...,[iV i /2] . (3.4) 

where [x] indicates the integer part of x. As mentioned above, V corresponds to the total 
number of quanta 

3 

V = $> , (3-5) 

i=i 

which is conserved by the interactions in (2.14) 
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The contributions to the Hamiltonian (2.14) involving 50(2) operators are diagonal 
in the basis (3.1) 

< [iVi], [N 2 ], [N 3 ];v 1 ,v 2 ,v 3 ;V\H l M \[N 1 ], [N 2 ], [N 3 }; Vl ,v 2 ,v 3 ;V >= -vj + N iVi ,(3.6) 

< [N 1 ],[N 2 ],[N 3 ] ] v 1 ,V2,v 3 ]V\H ij \[N 1 ],[N 2 ],N 3 y,v 1 ,V2,v 3 -,V > 

= 2v i v j -(v i N j + v j N i ) , (3.7) 

while the Vij operator has only non-diagonal matrix elements, since it involves the raising 
and lowering operators in the i,j indices, 

< [N,], [N 2 ], [N^v'^v^VlVijim, [N 2 ], [N 3 };v u v 2 ,v 3 ;V> 

= \jvj(v t + l)(Ni - v l )(N J - Vj + sv'.-.vi-i 

+ \Jvi(vj + l)(Nj - Vj)(Ni -Vi + l)5 v > ijV ._ 1 sv'j.vj+i > ( 3 - 8 ) 

with i,j = 1, 2, 3 and i < j. Note that because of the symmetry of the V 3 system the same 
number of bosons iVj = N (i.e. the same potential depth) is assigned to all three bonds. 
In the next Section the physical meaning of these interactions will become clear. These 
analytical results for the matrix elements and analogous ones for higher order interactions 
constitute one of the main advantages of the model. We point out, however, that although 
the local basis is convenient from a numerical point of view, it does not span the irreps of 
S 3 . A better way to carry out the diagonalization of (2.14) is to symmetrize the local basis 
(3.1), for which we can either symmetry- project the wave functions arising from the local 
basis once the Hamiltonian has been diagonalized, or generate the symmetry adapted 
one-phonon states and then construct the higher-phonon states by means of coupling 
coefficients. 13 '' For our purposes it is better to follow the second route, since in this way 
the wave functions explicitly carry the normal labels from the outset. In the case where 
spurious modes are present, the building-up procedure is essential, since in this way the 
unphysical modes can be exactly eliminated from the space. 
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In order to construct the normal basis we start by establishing explicit forms for the 
irreps of the group S3. For practical reasons it is convenient to work with real represen- 
tations, so we have chosen the cartesian harmonics as a basis for the E representation. 
In Table I we show the character table of the S3 group, including at the right the basis 
functions spanning the irreducible representations. In Table II we indicate the explicit 
irrep E carried by these functions in the reference frame of Figure 1. 

We now consider the one-phonon local functions. In this case the basis (3.1) has the 

form 

|l>=|[iV],[iV],[iV];100;l> , 
|2>=|[iV],[iV],[iV];010;l> , 

|3>=|[iV],[iV],[iV];001;l> , (3.12) 

which can be readily projected to the normalized states 

1 

Vf =^{2|l>-|2>-|3> 

= -Lj|2>-|3>} , (3.13) 

where we have used the notation for the wave functions. For higher phonon number 
the states are obtained through the coupling 13 ) 



V Al = A=< |1 > +|2 > +|3 > } , 



= E C( r i r 2r; 71727) Vl ^l V2 ^l , (3-14) 



Vl+V 2 

^7 

71,72 



where the coupled wave functions now correspond to a total number of phonons V = 
V1 + V2. The coupling (Clebsch-Gordan) coefficients C(; ) can be found in tables, 20 - ) or, in 
order to avoid phase inconsistencies, computed in a strightforward way using the explicit 
irrep given in Table II. 21 ) In Table III we present the Clebsch-Gordan coefficients derived 
in this fashion. 
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Using (3.14) repeatedly leads to a building-up procedure to derive the symmetry 
adapted basis for higher-phonon numbers. To achieve this task, however, we must obtain 
the decomposition of the products (Ai) VA i (E) VE , where vau ve correspond to the number 
of phonons in the normal modes. In Table IV we indicate, as an example, the reductions 
for two and three quanta. The procedure to obtain these reductions is a standard one, 
explained in many group theory textbooks. 21 ) One should bear in mind, however, that 
the boson nature of the vibrations implies that only symmetrized products are allowed 
for phonons in the same mode. For example, for two quanta in the E mode we have 
in general the following reduction in terms of the symmetric [ ] and antisymmetric { } 
contributions 21 ) 

E <g> E = [E <g> E]@ {E <g> E} , 

where 

[E ® E] = A 1 © E ; {E <g> E} = A 2 . 

Since the two phonon state (3.14) associated to the product {E®E} vanishes automatically, 
we are left only with the symmetrized product [E ® E]. 

The general procedure is now clear. Once the form of the Hamiltonian has been deter- 
mined by symmetry considerations, we proceed to construct the symmetry adapted basis 
by projecting the one-phonon local functions. The higher-phonon functions are then gen- 
erated from the one-phonon symmetrized states by means of the coupling (3.14). Finally, 
we carry out the diagonalization in the symmetrized basis, where full advantage can be 
taken of group-theoretical properties. In particular, the Hamiltonian matrix separates into 
blocks corresponding to the irreps of the symmetry group S3. For example, from Table 
IV we see that in the three phonon manifold the number of functions is 10, which reduce 
to three blocks of dimensions 3x3, lxl and 3x3, corresponding to the irreps A±, A2 
and E, respectively. The simplification becomes more significant as the complexity of the 
molecular system and/or the phonon number increase. 13 ) 
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4. Analysis of Interactions 



We now proceed to analyze the interactions involved in the Hamiltonian (2.14). In 
order to do so it is convenient to recall the standard labeling of states as well as the 
vibrational Dunham expansion for T^-triatomic molecules. 2 ) 

As is well known, this type of molecules exhibit three vibrational degrees of freedom, 
which give rise to two normal modes associated to A\ and E symmetries. The normal 
states are then specified by the number of quanta in each mode \v Al , ve >• In addition, 
the double degenerate E mode carries an intrinsic angular momentum /, whose values 
depend on % and are given by 2 - 1 

l = v E , ve-2,...,1 or . (4.1) 

The states are then specified by the quatum numbers v A , ve and / with the notation 

\v Al , v l E > . (4.2) 

The simplest way to reproduce the general features of the spectrum is by means of a 
Dunham expansion, which up to quadratic terms takes the form 2 ) 

E v (v Al ,v E ,l) =E + u^(v Al + l/2) + uf(v E + l) 

- XeUjf 1 ( VAl + 1/2) 2 - X e U)f{v E + I)' 

+ x 12 (v A + l/2)(v E + l) +922I 2 , (4.3) 

The first two terms in the sum correspond to the harmonic contributions to the energy, 
while the next three terms, (v A + 1/2) 2 , (ve + l) 2 and (v Al + 1/2) (ve + 1), represent 
the first anharmonic corrections. The last term is the intrinsic (or vibrational) angular 
momentum and gives rise to the correct ordering for states with the same value of the 
quantum numbers v A and ve- In addition, each wave function \v Al , v l E > carries a definite 
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symmetry, which is closely related to the / quantum number. For / = the states are 
totally symmetric and labeled as a\ (we use lower case letters to denote the symmetry). 
For I = 3/c, k = 1, 2, 3, . . . , two levels corresponding to a\ and ai symmetries appear, while 
for I = 3k + 1 or I = 3k + 2, k = 0, 1, 2, . . ., the states exhibit e symmetry. Note that 
for / 7^ there are two components ±/, although this is not explicit in the notation. 2 ) We 
remark that the expansion (4.3) does not remove the degeneracy of the levels a\ and ci2 
associated to the / = 3k states. The same is true for any order in the Dunham expansion. 
Experimentally this degeneracy is not present, but this cannot be taken into account by 
such simple parametrizations. 

Let us now analyze the interactions involved in (2.14). If we compute the matrix 
elements of the operators ^ij an d Yl Vij m the one-phonon manifold, we obtain 



K 1 V^lX^J'lV 41 >= ~ 2N ' ( 4 - 4a ) 

< 1 ^ Al \^V ij \ 1 ^ Al >=2N , (4.46) 

K 1 VflX^/Vf >=-2iV , 7 = 1,2 . (4.5a) 

< 1 <|^^ J | 1 <>=-iV , 7 = 1,2. (4.56) 

From these results we conclude that the operator 

W B = _I{£V y + X;fry} (4.6) 



does not contribute to the energy of the A\ mode. The —1/3 factor was added for later 
convenience. In Figure 2 we show the spectrum generated by (4.6) for the two and three 
phonon manifolds, as a function of N. Note that for large N the operator (4.6) behaves 
as fiEi the number of phonons in the E mode, a result we shall explain in Section 5. 

The previous analysis leads to the question of whether it is possible to construct 
from (2.15) an operator affecting only the A\ mode. This is indeed possible and through 
projection we find 

WAx^^Vy-^fly} , (4-7) 
15 



which satisfies 

< 1 ^f\n Al \ 1 ^>=0 , 7 = 1,2 , (4.8) 

as required. We show in Figure 3 the spectrum generated by (4.7) as a function of N. For 
large N the operator (4.7) behaves as fiA i: the number of phonons in the A\ mode (see 
Section 5). 

We have thus constructed operators that selectively affect to the A\ and E modes. 
The Hamiltonian (2.4), however, includes three independent operators. A third operator 
can be easily derived: 

i=l i,j = l 

which is diagonal in the local basis (3.1) and satisfies 

< 1 >= , (4.10a) 

< 1 ^\V\ 1 ^>=0, 7 = 1,2. (4.106) 

In Figure 4 we schematically show the effect of V as a function of N in the two and three- 
phonon manifolds. Note that this operator is diagonal in the local basis and vanishes in 
the large N limit. 

The method followed in this section is general and indicates a procedure to define 
operators with definite actions over the physical space. Additionally, the use of these 
symmetry adapted operators significantly improves the convergence of the mean square 
search of parameters in the diagonalization of the Hamiltonian. 

Although the operators (4.6), (4.7) and (4.9) induce the characteristics of the spectrum 
generated by the harmonic and anharmonic contributions in (4.3), they cannot reproduce 
the effect of the Z 2 term. This term orders, for a positive (negative) value of (722, the 
vibrational levels in each phonon multiplet (va 1 ,v 1 e ) according to increasing (decreasing) 
value of /. Figures 2-4 indicate that the characteristic pattern of the V^h vibrational 
spectrum cannot be reproduced by the simple Hamiltonian 

H = aH Al +(3H E + lV . (4.11) 
16 



This analysis, however, does show that the simple algebraic Hamiltonian (2.14) can be 
interpreted in a physically meaningful way by concentrating on the symmetry properties 
of the interactions, as expressed in (4.11). In the next section we present a systematic 
procedure to derive the full set of interactions in the algebraic framework, starting from 
those present in configuration space. 



5. Algebraic Interactions and Configuration Space Operators 

In order to establish the algebraic representation of configuration-space operators, we 
start by analyzing the harmonic limit of the angular momentum operators 

[J , J±] =±4 , (5.1a) 
[J+,J_]=2J . (5.16) 

The action of the J± on the Morse states | [N] , v > is given by 

J+\[N],v>= a/ v(N — v + 1) \[N],v-l> , (5.2a) 



J-\[N],v>=y/(N-v)(v + l)\[N],v + l> , (5.26) 
where N and v were defined in (3.2) and (3.3). Defining the change of scale transformation 

b=^= , 6t = ^= , (5.3) 

it is clear that 



lim b \[N]v >= y/v\[N],v- 1 > , (5.4a) 
lim tf\[N]v >= + 1 \[N],v+ 1 > . (5.46) 

which correspond to the harmonic limit of the model, as expected from the role of iV in 
Eq. (2.4), i.e., for infinite potential depth the Morse potential cannot be distinguished 
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from an harmonic potential. Using the definitions (3.2), (3.3) and (5.3), we can rewrite 
the commutation relation (5.1a) in the new form: 

[S,6t] = l-^ , (5.5a) 



where 



N 

v=--J (5.56) 



is the Morse phonon operator corresponding to the definition (3.2). The limit N — > oo 
leads to the contraction of the SU(2) algebra to the Weyl algebra generated by 6, and 
1, with the usual boson commutation relation [6,6'] = 1. Relations (5.3-5.5) indicate 
the procedure to arrive at the harmonic limit of the model. Each J +i , J_i should be 
renormalized by dividing by ^/Wl and then take the limit iVj — > oo. 

As an example of this procedure we take the harmonic limit of the Morse Hamiltonian 

(2.7) 



lim —H M = lim — 

N^ooN N^oo N 



\ (j_J + + J + J_) - 



N 
~2 



= I( 6 t6 + 66t)-i 

= 6 f 6 , (5.6) 

which has eigenvalues rib, i n agreement with the harmonic limit of Eq. (2.9). Applying 
the same procedure to the symmetry projected interactions of Eq. (4.11) we find 

lim — D = , (5.7) 

where fiA 1 and n e are the operators corresponding to the number of phonons in the A\ and 
the E modes, respectively, as can be readly shown using the technique discussed below. 
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We can now interpret Eq. (5.3) in the opposite sense, i.e. as a way to construct 
the anharmonic representation of harmonic operators. Any given function of 6, &t can be 
mapped into the same function of J+, J_ through the correspondence 

^ % , 6 f - "4 • (5-8) 



As an example, we consider again the one-dimensional harmonic oscillator 

H=^(tfb + btf) , (5.9) 

with eigenvalues E = rib + 1/2, and follow the reverse order in (5.6). To obtain its 
anharmonic representation we carry out the correspondence (5.8) to get 

= ±(j 2 -Ji)=v + l/2-v 2 /N , (5.10) 

which is the algebraic realization of the Morse oscillator, as shown in Section 2. 

The general procedure to derive the algebraic realization of a given configuration- 
space operator is thus the following. We first write down the operator in terms of normal 
coordinates and momenta {g, p}^, and express it in terms of the harmonic bosons 

b^ = ^=(q + ip)^ , (5.11) 

where V denotes the irrep spanned by the normal coordinate and 7 is its row. In the next 
step we write down the normal operators {6p^,&r 7 } m terms of local ones {b\ 7 bi}. The 
explicit relations are of the general form 



i 

b r 7 =J2^,A , (5.12) 



1 
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where the set {ct^} can be obtained by projecting the local operators {b\} on the irrep 
spanned by the normal coordinates {q,p}^. Finally, we substitute (5.12) in the expression 
of the interaction given in terms of the normal operators (5.11) and carry out the corre- 
spondence (5.8). As an example of this procedure we derive the algebraic form of the Z 2 
interaction in (4.3). 

The representation of the operator I in terms of normal coordinates qf , q 2 , which 
span the irrep E, is given by 22 ) 



e d e d 



''^Bg'^m ■ (B13) 

Introducing the operators (5.11), this expression transforms into 

l = -i(b?b?-b?h E ) . (5.14) 

We now write the normal operators {bf\b^} in terms of the local ones {b\}. This can 
be done using the projected functions (3.13) 

b? = ±={2b\-b\-b\) , 

^ = ^(4-4) , (5-15) 

and equivalent expressions for the annihilation operators. Finally, we substitute (5.15) 
into (5.14) to obtain 

i A2 = -^{b\(b 2 -h)+bUb 3 -b 1 )+bUb 1 -b 2 )} , (5.16) 

where we have explicitly indicated the irrep carried by the / operator. The fact that it 
corresponds to an A 2 symmetry can be deduced either by analyzing the transformation of 
I under the £3 group or by identifying the Clebsch-Gordan coefficients C(EEA 2 ] 7172I) 
in (5.14). The corresponding realization in the model is then obtained by applying the 
correspondence (5.8) 

U 2 = -^{J.^j+a - J +3 ) + J_ 2 (J +3 - + J_ 3 (J +1 - J +2 )} . (5.17) 
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In turn, the l\ operator is obtained by squaring (5.17). The same kind of analysis can be 
applied to arbitrary configuration space interactions. 22 ) 

We have presented in this section a general method to derive the realization of op- 
erators in the algebraic model, starting from their representation in configuration space. 
This procedure considerably increases the power of the algebraic approach, since it can 
be used to incorporate into the model the fundamental interactions known from the con- 
figuration space methods. Note that this procedure allows, in principle, to establish the 
relation between the algebraic parameters and the force field strengths obtained from ab 
initio calculations. It is also possible, however, to apply the model in a purely algebraic 
fashion and still deduce the fundamental interactions, as we explain in the next Section. 



6. The Anharmonic Oscillator Symmetry Model 

In this Section we present a general framework to construct all interactions in the alge- 
braic model in a systematic way. We shall henceforth refer to this procedure, together with 
the methods introduced in the previous sections, as the Anharmonic Oscillator Symmetry 
Model (AOSM). 

We start by introducing a set of generators with well-defined tensorial properties under 
the point group 

A»,7 = a 7,i AM ' (6-1) 
i 

where \i = +, — , 0. For the case of T>zh molecules the expansion coefficients are the same 
as those in Eq. (5.12). We then construct from these symmetry projected generators a set 
of interactions that are scalars under the point group, such as 

T t (T) = \ £ {J T -Jl n + JljQ , (6.2a) 

7 

and 
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t Ai = (r) J2 j oJo,t ■ ( 6 - 2& ) 

7 

Higher order tensors can be systematically constructed by means of (6.1) and the Clebsch- 
Gordan coefficients for the point group. 

For triatomic I^-molecules V = Ai,E, the relevant symmetry projected generators 

are 



}e _ 



V3 
1 

71 



fE . 



J= (^,2 - 4,3) , (6.3) 



with fj, = +, — ,0. According to (6.2) we can construct four possible interactions that are 
quadratic in «/£ 7 , three of which correspond to linear combinations of the terms in (2.14) 

3 „ 3 



3 3 



(Ai) = ?iV 2 -(£H™- £ #y) , (6.4c) 

«=i »>j=i 

while the fourth is not independent, since 

3 
i=l 

In addition to the operators in (6.2) which transform as A\ under D3, we can also construct 
other bilinear combinations with well-defined tensor properties, 

rpE _ 1 ( fE fE fE fE \ 

rpE _ 1 / ?E 7 E , ?E tE \ 
J 2 - g l v J -,l J + ,2 + J -,2 J +,lj ' 

^^(i-V + £ 2^Vf,i) • (6-5) 
22 



The operator T A2 is proportional to the intrinsic angular momentum operator of (5.17). 

In lowest order (i.e. quadratic in Tf ) there are two possible A\ interactions, (Tf) 2 + 
(Tf ) 2 and (T A2 ) 2 . In order to interpret these interactions we take the harmonic limit of 
(6.5) 

Jim ^ = L z = l - (tf E b E2 - tf E b El ) . (6.6) 

The operators L x , L y and L z close under the commutation relations of SU(2). In the 
harmonic limit (T^ 2 ) 2 corresponds to L 2 = Z 2 /4 (see eq. (5.14)) while (Tf) 2 + (Tf ) 2 
goes to L 2 + L 2 = L 2 — L 2 with I? = tie^e + 2)/4, corresponding to an anharmonic 
contribution to the E mode. 

In the next order (i.e. cubic in T^) we can first couple Tf to E and then couple again 
the resulting operator to Tf to obtain an A\ operator, 

6 M = (ff Tf - ffff^j Tf + (if + ff Tf ) Tf 

= ~(Tj + Tf) , (6.7) 

where we have introduced T± = Tf ± zT 2 E , which in the harmonic limit reduce to the 
ladder operators L± = L x ± iL y . In the harmonic limit the operator (6.8) has a clear 
physical interpretation: it couples states with AL Z = ±3, or expressed in terms of the 
intrinsic angular momentum, it couples states with the same v E and Al = ±6. For this 
reason it splits the ai and vibrations that are associated with the (va 1 , v 1 ^ 3 ) multiplet. 

The explicit realization of Al in the AOSM can be obtained by expressing T + (and 
T_ = (T + )t) in terms of the J± t i operators through equations (6.5) and (6.3). The final 
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result is 



T + = 4 



6 



+ 



2J_i J_|_i — J_ 2 J+2 — J-3J+3 

— J-l(J+2 + J+3) — J-2{J+l ~ 2J+3 

— J-s(J+i — 2J+2) 

J-l(J+2 — J+3) + J-2(J+1 — J+2) 



2V3 _ 

+ J_ 3 (J+3 - J+l) 



(6.8) 



In the next section we show that the interaction (6.7) is essential to describe the highly 
anharmonic molecule . It should be clear that the AOSM can be applied in a similar 
way to molecules exhibiting arbitrary symmetry groups. We have thus presented a sys- 
tematic procedure to construct, up to a certain order, all relevant interactions, based on 
the introduction of operators with well-defined tensorial properties under the point group 
(see e.g. equation (6.1)), which can then be combined into Vsh scalar interactions. 



7. Application to HJ, Be 3 and Na^ 

In this section we apply the AOSM to Na^~, Be 3 and H^j". As mentioned in the 
Introduction, we have chosen these molecules because they exhibit a wide range of behavior, 
ranging from the very anharmonic spectrum of to the almost exact harmonic spectrum 
of Na^. 

According to the discussion presented in the previous section, the Hamiltonian 

H = aH Al + (3H E + 7V + 51 2 , (7.1) 

contains the main physical interactions that describe a D3/,-triatomic molecule, whose 
spectrum is close to the one generated by the Dunham expansion (4.3). As mentioned 
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before, the Dunham expansion implies a degeneracy between the ai and a<i levels associated 
to the same quantum number /, while the spectrum generated by (7.1) does lead to their 
splitting, although it is generally small. Experimentally this splitting is observed, even for 



effect is produced by the 7Yai operator, as can be seen from Fig. 3 . 

We now consider the molecules Hjj~, Be3 and NaJ. The atoms in the first molecule 
are very light and the spectrum is highly anharmonic, a fact that is reflected by a strong 
splitting of the ai, a2 levels (200 cm -1 for V = 3), as well as by a relatively large splitting 
of levels belonging to the same (va 1 , ve) multiplet but with different value of / (220 cm -1 
for V = 2). This is in contrast with the case of Na^, where the splitting between the 
ai and levels is very small (0.11 cm -1 for V = 3) as well as the /-dependent splitting 
(0.82 cm -1 for V = 2). The molecule Be3 exhibits an intermediate behavior, although in 
this case an a\, splitting is not present in the fitted data, since we have generated its 
spectrum from an ah initio calculation where no splitting terms are included. 15 '' 

In Table V we present a least square fit calculation for H^j~, Be3 and Najj" up to three 
quanta, using the Hamiltonian (7.1) with 5 = 0. The standard deviation (rms) was taken 
to be 



where n and n p correspond to the number of fitted levels and parameters involved, respec- 
tively. From this calculation we find a large difference in the quality of the fit between Hg~ 
and the other two molecules. In Table VI we present the same calculations, but including 
the I 2 interaction. We see that the difference in quality persists; while the Hamiltonian 



of Carter and Meyer, 16 ) who are forced to include twice as many terms in the potential 
for H^j~ than for the Na^ molecule. The simplest possible set of such interactions in the 



molecules like Nag", in spite of its harmonic behavior. In the model Hamiltonian (7.1) this 




(7.2) 



(7.1) is quite sufficient to describe the Be3 and NaJ molecules, we clearly require ad- 
ditional interactions to properly describe H^ . This fact is in accordance with the work 



25 



AOSM correspond to higher powers of the symmetry adapted operators (4.6) and (4.7). 
We propose 

n E> rtA x E = ^ • \ L6 > 

If we add this set of interactions to the Hamiltonian (7.1) in the energy fit for H+ the 
rms deviation reduces to 15.74 cm -1 . It is possible to further improve the fit by taking 
into account the interaction (6.7) in addition to the set (7.3). A more general algebraic 
Hamiltonian to describe V^h molecules is then 

H = aH Al + PHe + 7V + 5P + a [2] H 2 Ai 

+ 0n 2 E +^n Al E+e(fl+T^) . (7.4) 

As mentioned before, the operator T+ + T£ has the effect of splitting the a\ and ai levels 
arising from the same angular momentum /, which explains the need for this interaction in 
Hg". In Table VII we present the least square energy fit to Hg" using the Hamiltonian (7.4), 
with an rms deviation of 5.84 cm -1 . We remark that in order to describe H^ for higher 
phonon numbers we need to include higher order interactions. We believe, however, that 
this result is very encouraging. If we omit the purely anharmonic interaction V the rms 
increases to 24.37 cm -1 , while carrying out the harmonic limit (N — > oo), where V = 0, 
the rms obtained is 31.17 cm -1 . 



8. Conclusions 

We have introduced the Anharmonic Oscillator Symmetry Model and applied it to a 
set of Dsft-triatomic molecules. The model is based on symmetry methods which systemat- 
ically incorporate group theoretical techniques, providing a clear methodological procedure 
that can be applied to more complex molecules. We have introduced symmetry adapted 
operators that have a specific action over the function space. This is a general procedure 
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which gives rise to a clear physical interpretation of the interactions and has the addi- 
tional advantage of considerably improving the convergence during the least square energy 
fit. Furthermore, based on the harmonic limit of the SU(2) algebra we have proposed a 
systematic approach to derive an algebraic realization of interactions given in configura- 
tion space. The model surmounts one of the main objections raised against the use of 
algebraic models, where it was not possible to obtain a direct correspondence with the 
configuration-space approaches. Although we have illustrated this procedure by means of 
the I 2 interaction, the method can be used for arbitrary operators. For the general case 
when there is no information about the form of these interactions in configuration space, 
we have devised an algebraic procedure to derive them using their tensorial structure un- 
der the point group. The combination of the different methodologies leads to the AOSM, 
which can be applied in the same fashion to more complex molecules. We remark that 
the model can be extended in several ways. For example, Fermi resonances can be taking 
into account using perturbation theory, while the rotational degrees of freedom can be 
incorporated by coupling the vibrational wave functions to rotational states carrying the 
appropriate point symmetries. 23 ' 24 ) 

We believe that the AOSM represents a systematic, simple but accurate alternative 
to configuration space methods, particularly for polyatomic molecules, where the integro- 
differential approaches are too complex to be applied or require very large numerical cal- 
culations. Since the model provides manageable wave functions, it is possible to evaluate 
the matrix elements of arbitrary physical operators, which have a simple representation in 
the algebraic space. A finer test for the model is to use these wave functions, for exam- 
ple, for the evaluation of infrared and Raman intensities. The transition operators can be 
constructed by applying our method to the configuration-space parametrizations, which 
correspond to the harmonic limit, N — > oo, or purely algebraically by using their tensorial 
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properties under the corresponding point group. The analysis of electromagnetic intensi- 
ties, as well as the application of the model to other molecular systems will be presented 
in future publications. 17 ' 23 ) 
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Figure Captions 



Figure 1. Assignment of the U l (2) algebras for V^h triatomic molecules and the selection of 
the Cartesian coordinate system. The elements of the symmetry group V3 are also 
indicated. 

Figure 2. Eigenvalues of the operator H E as a function of the number of bosons N for (a) two 
phonons and (b) three phonons. 

Figure 3. Eigenvalues of the operator as a function of the number of bosons N for (a) two 
phonons and (b) three phonons. 

Figure 4. Eigenvalues of the operator V as a function of the number of bosons N for two and 
three phonons. 
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Table I. Character table for the V 3 pa S3 group. The set {R x , R y , R z } represents the 
components of an axial vector. 



V 3 


E 


C 3 (2) C 2 (3) 


Basis functions 


A x 


1 


1 1 


— rj n rj 

2z z — x z — y z 


A 2 


1 


1 -1 


z- R 


E 


2 


-1 


(x,y) , (R x ,R y ); 








(zx,yz) (xy,x 2 -y 2 ) 



Table II. The E irreducible representation of the generators of the £3 group in the 

Cartesian basis of Table I. 

Irrep (123) (23) S 3 
C3 V3 



E 




Table III. Clebsch-Gordan coefficients C(fjivT;ij'y) for the £3 group compatible with the 
irreps of Table II. 

A 2 xE (11) (12) 



(E) 1 1 



(r) 7 C{nvY;iji) (£;)2 _ 1 Q 



E x E 


(11) 


(12) 


(21) 


(22) 


Ax 


1/V2 








1/V2 


A 2 





1/V2 


-1/V2 





(E)i 


1/V2 








-1/V2 


(Eh 





-1/V2 


-1/V2 
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Table IV. V 3 decompositions for two and three quanta. 

A\ E Dimension of the 

va 1 ve Irreps reducible representation 



2 





A 1 




1 





2 




BE 


3 


1 


1 


E 




2 










6 


3 





A x 




1 





3 


Ati 


D A 2 ®E 


4 


2 


1 


E 




2 


1 


2 




BE 


3 



10 



32 



Table V. Least square energy fit for H3", Be 3 and Na^. All energies are in cm -1 . We 
indicate the rms deviation (7.2) and the parameters obtained. The number of bosons N 
was taken to be 30. The energy difference is given by AE = E t h — E exp , where E exp are 
taken from ab initio calculations. 



(va^e) Symmetry 



(01 1 ) 


e 


(10°) 


ai 


(02°) 
(02 2 ) 


ai 


e 


(H 1 ) 


e 


(20°) 


ai 


(03 1 ) 


e 


(03 3 ) 


ai 


(03 3 ) 


a 2 


(12°) 


ai 


(12 2 ) 


e 


(21 1 ) 


e 


(30°) 


ai 



Exp. 16) AE 

2521.27 -0.17 

3178.32 -29.32 

4777.02 102.26 

4997.41 -45.57 
5553.67 -14.97 
6261.92 -22.77 

7003.49 105.59 

7282.52 -102.30 

7492.64 -105.84 

7769.09 116.14 

7868.64 -32.58 

8486.90 1.37 

9251.42 14.89 



Be 3 

Exp. 15 ') AE 

399.10 4.09 

458.40 0.25 

782.40 8.33 

794.40 1.31 

845.10 2.02 

907.60 -0.04 

1161.90 3.91 
1185.90 -13.85 

1185.90 -3.95 

1216.00 9.71 

1228.00 -3.37 

1281.90 -0.72 

1347.60 -1.04 



Na+ 

Exp. 16 ) AE 

99.95 1.68 

140.45 2.16 

198.90 2.02 

199.72 1.25 

239.29 1.70 

280.35 1.66 

297.67 0.25 

299.26 -1.25 

297.67 -1.36 

337.19 -0.15 

337.94 -0.86 

378.06 -0.91 

419.70 -1.49 



rms 78.55 6.48 2.58 

a 3148.996 458.653 142.608 

Parameters (3 2521.105 403.185 101.633 

7 3796.387 328.843 49.434 
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Table VI. Least square energy fit for H+, Be 3 and NaJ using the Hamiltonian (7.1). We 
show the energy differences AE = E th — E exp . The values of the energies E exp are given 
in Table V. 









Be 3 


Na+ 


/ / \ 

{VAy E ) 


Symmetry 


AE 


AE 


AE 


(01 ) 


e 


-37.18 


0.51 


0.93 


(10°) 


at 


-21.70 


0.02 


1.95 


(02°) 


a x 


-16.32 


-0.74 


0.37 


(02 2 ) 


e 


-33.84 


0.17 


0.84 


(H 1 ) 


e 


-35.74 


0.82 


1.68 


(20°) 


ax 


-13.42 


-0.04 


1.26 


(03 ) 


e 


18.66 


-2.05 


—1.19 


(03 ) 


a x 


16.05 


— 1.23 


-0.34 


(03 ) 


a 2 


-0.62 


0.61 


-0.33 


(12°) 


ax 


46.62 


1.90 


-0.01 


(12 2 ) 


e 


-12.83 


-1.36 


0.34 


(21 1 ) 


e 


-3.38 


0.79 


-0.19 


(30°) 


ax 


22.56 


-1.66 


-2.06 




rms 


30.15 


1.24 


1.25 




a 


3156.616 


458.911 


142.396 


Parameters 




2446.638 


396.265 


100.317 




7 


3131.825 


209.744 


21.312 




S 


-12.485 


-0.9533 


-0.1867 
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Table VII. Least square energy fit of the H3" molecule using the Hamiltonian (7.6). 



(va 1 v l E ) Symmetry AE 



(01 1 ) 


e 


-1.55 


(10°) 


ai 


0.42 


(02°) 


ai 


7.48 


(02 2 ) 


e 


-5.69 


( n ) 


e 


n a 1 
— O.Dl 


(20°) 


ai 


-0.11 


(03 1 ) 


e 


-4.46 


(03 3 ) 


ai 


3.18 


(03 3 ) 


a 2 


2.44 


(12°) 


ai 


0.66 


(12 2 ) 


e 


-5.0 


(21 1 ) 


e 


4.07 


(30°) 


ai 


-1.23 


rms 




5.84 cm' 1 



Parameters 



a = 3193.600 
(3 = 2507.157 
7 = 2807.833 



5 = -13.439 

a® = -14.855 
pm = -27.752 
f [2] = -28.043 



e = -0.900 
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